Ideal strengths and bonding properties of PuC>2 under tension 
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Abstract 

We perform a first-principles computational tensile test on PuC>2 based on density-functional 
theory within local density approximation (LDA)+£/ formalism to investigate its structural, me- 
chanical, magnetic, and intrinsic bonding properties in the four representative directions: [001], 
[100], [110], and [111]. The stress-strain relations show that the ideal tensile strengths in the four 
directions are 81.2, 80.5, 28.3, and 16.8 GPa at strains of 0.36, 0.36, 0.22, and 0.18, respectively. 
The [001] and [100] directions are prominently stronger than other two directions since that more 
Pu— O bonds participate in the pulling process. Through charge and density of states analysis 
along the [001] direction, we find that the strong mixed ionic/covalent character of Pu— O bond 
is weakened by tensile strain and Pu02 will exhibit an insulator-to-metal transition after tensile 
stress exceeds about 79 GPa. 
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I. INTRODUCTION 



Plutonium-based materials have been extensively studied due to their interesting physical 
behaviors of the 5/ states and have always attracted particular attention because of their 
importance in nuclear fuel cycle [l-4]. Recently, experimental reports 5-7] on the strategies 
of storage of Pu-based waste illustrated that metallic plutonium surface easily oxidizes to 
PU2O3 and PuC>2 in surrounding of air and moisture. Previous calculations Q using the hy- 
brid density functional of (Heyd, Scuseria, and Enzerhof ) HSE have established trends of the 
electronic properties of actinide dioxides AnC>2 (An=Th— Es). Along the series, Mott insula- 
tors (/—>•/) are obtained prior to PuC>2, whereas above AH1O2, charge-transfer insulators (O 
2p — ^An 5/) are observed. PUO2 and Am02 lie in a crossover from localized to delocalized 
/ electron character with increasing Z. They have strong 5f-2p orbital energy degeneracy, 
which can lead to unexpected orbital mixing. Recent theoretical work performed by Petit 
et al [8j demonstrated that the dioxide is the most stable oxide for the actinides from Np 
onward. All these works indicate that investigations of PuC>2 have particular meaning in 
actinide compounds. 

At ambient condition, PuC>2 crystallizes in a fluorite structure with space group Fm3>m. 
And at 39 GPa, PuC>2 undergoes a phase transition to an orthorhombic structure of cotun- 



nite type with space group Pnma [9]. In our previous work [1CH12|. we have systematically 
investigated the structural, electronic, mechanical, and optical properties of fluorite PuC>2 
within the LDA+ U and GGA+ U formalisms. We found that comparing with the experi- 
mental data and the theoretical results, the accuracy of our atomic-structure prediction for 
antiferromagnetic (AFM) PuC>2 is quite satisfactory by tuning the effective Hubbard param- 
eter U in a range of 3-4 eV within the LDA/GGA+f/ approaches. Subsequently, Jomard 
et al [3] confirmed our results and further reported optical and thermodynamic properties 
of Pu02- In 2008, Yin and Savrasov [l4[ successfully obtained the phonon dispersions of 
both UO2 and Pu02 by employing the LDA+DMFT scheme. In 2009, Minamoto et al. 
is) ] investigated the thermodynamic properties of Pu02 based on their calculated phonon 
dispersion within the pure GGA scheme. 

However, although many efforts have been performed on PUO2, little is known on its 



theoretical tensile strength 



llj . The ideal strength of materials is the stress that is re- 



quired to force deformation or fracture at the elastic instability. Although the strength of 
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a real crystal can be changed by the existing cracks, dislocations, grain boundaries, and 
other microstructural features, its theoretical value can never be raised, i.e., the theoretical 
strength sets an upper bound on the attainable stress. In this study, we focus our sight on 
the structural, electronic, mechanical, and magnetic features of PuC>2 under tension. The 
stress-strain relationships and the ideal tensile strengths are obtained by performing a first- 



principles computational tensile test (FPCTT) 16] . Under tension, the bonding nature and 



electronic occupation characters are systematically studied. 



II. COMPUTATIONAL METHODS 



The first-principles density-functional theory (DFT) ca 



culations on the basis of the 



frozen-core projected augmented wave method of Blochl [17] are performed within the Vi- 
enna ab initio simulation package (VASP) [181 ]. where the LDA 19) is employed to describe 
electron exchange and correlation. For the plane- wave set, a cutoff energy of 500 eV is 
used. The &-point mesh in the full wed ge o f the Brillouin zone (BZ) is sampled by 9x9x9 



grids according to the Monkhorst-Pack [20j for the fluorite unit cell and all atoms are fully 
relaxed until the Hellmann-Feynman forces become less than 0.001 eV/A. The plutonium 
Qs 2 7s 2 6p 6 6d 2 5f A and the oxygen 2s 2 2p 4 electrons are treated as valence electrons. The 
strong on-site Coulomb repulsion among the localized Pu 5/ electrons is described by using 



the LDA+ U formalism formulated by Dudarev et al. [2l{| . In this paper the Coulomb en- 
ergy (U) and the exchange energy (J) are set to be constants: [7=4.75 eV and J=0.75 eV. 
These values of U and J are the same as those in our previous study of plutonium oxides 



10 



11] . Using these parameters, the LDA+U gives ao=5.362 A, which is very close to the 
experimental value of 5.398 A [5j. And our results reproduce all the features included in 
our previous work [lOj. In particular, we recover the main conclusion that although the 
pure LDA fail to depict the electronic structure, especially the insulating nature and the 
occupied-state character of Pu02, by tuning the effective Hubbard parameter in a reasonable 
range, the LDA+ U approaches can prominently improve upon the pure LDA calculations 
and, thus, can provide a satisfactory qualitative electronic structure description comparable 
with the photoemission experiments f], Spin-polarized calculations are performed and 
we find that the AFM spin alignment is the most stable configuration among nonmagnetic, 
ferromagnetic (FM), and AFM configurations. The total-energy difference (E FM —E\ FM per 
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formula unit at respective optimum geometries) within the LDA+ U is calculated to be 0.705 
eV. 

In the FPCTT, the stress-strain relationship and the ideal tensile strength are calcu- 
lated by deforming PuC>2 crystal to failure. The anisotropy of the tensile strength is tested 
by pulling the initial ffuorite structure along the [001], [100], [110], and [111] directions. 
As shown in Fig. [TJ three geometric structures are constructed to investigate the tensile 
strengths in the four typical crystallographic directions: [T](a) shows a general ffuorite struc- 
ture of PuC>2 with four Pu and eight O atoms; [TJ[b) a body-centered tetrahedral (bet) unitcell 
with two Pu and four O; and [Tfc) a orthorhombic unitcell with six Pu and twelve O. In 
FPCTT, the tensile stress is calculated according to the Nielsen-Martin scheme [22] 

1 dE tota \ 

'« = n aiv (1) 

where E a $ is the strain tensor (a, (3=1,2,3) and Q is the volume at the given tensile strain. 
Tensile processes along the [001], [100], [110], and [111] directions are implemented by in- 
creasing the lattice constants of these three orientations step by step. At each step, the 
structure is fully relaxed until all other five stress components vanish except that in the 
tensile direction. 



III. RESULTS 

1. Theoretical tensile strength 

The calculated total energy, stress, and spin moments as functions of uniaxial tensile 
strain for AFM Pu0 2 in the [001], [100], [110], and [111] directions are shown in Fig. El 
Evolutions of the lattice parameters with strain in all four tensile processes are presented 
in Fig. [3j Clearly, all four energy-strain curves increase with increasing tensile strain, but 
one can easily find the inflexions by performing differentiations. Actually, at strains of 0.36, 
0.36, 0.22, and 0.18, the stresses reach maxima of 81.2, 80.5, 28.3, and 16.8 GPa for pulling 
in the [001], [100], [110], and [111] directions, respectively. These results clearly indicate 
that the [001] direction is the strongest tensile direction and [111] the weakest. 

In fact, there are eight Pu— O bonds per formula unit for fluorite Pu02- The angle of 
all eight bonds with respect to the pulling direction is 45° in [001] direction. For pulling 
direction of [100], the bonding structure is same with that of [001] direction. So their 
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FIG. 1: (Color online) Schematic illustration of tension along (a) [001] or [100], (b) [110], and (c) 
[111] orientations. Blue atoms are plutonium atoms while white atoms are oxygen atoms. The 
AFM order is indicated by white arrows. 
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FIG. 2: (Color online) Dependence of the (a) total energy (per formula unit), (b) stress, and (c) 
spin moments on tensile strain for AFM Pu02 in the [001], [100], [110], and [111] directions. 

pressure behaviors of total energy, stress, and spin moments are almost same before strain 
of 0.40. The special behaviors of total energy, stress, and spin moments over e=0.40 for 
pulling in the [100] direction have tight relations with the abrupt behavior of its crystal 
structure. In fact, over e=0.40 the face-centered tetragonal (fct) structure will transit into an 
orthorhombic structure, as indicated by Fig. [3]^b). The lattice parameter of [001] direction is 
shortened and [010] direction elongated by tensile deformation. However, in [110] direction 
only four bonds make an angle of 45° with the pulling direction. Four other bonds are 
vertical to the pulling direction. In [111] direction, two bonds are parallel to the pulling 
direction and six others make an angle of about 19.5° with the pulling direction. It is easy 
to understand that the bonds vertical to the pulling direction have no contributions on the 
tensile strength and the bonds parallel to the pulling direction are easy to fracture under 
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FIG. 3: (Color online) Dependence of the lattice parameters on tensile strain for AFM PuG"2 in 
the (a) [001], (b) [100], (c) [110], and (d) [111] directions. 

tensile deformation. Therefore, the fact that the tensile strength along the [001] or [001] 
direction is stronger than that along [110] and [111] directions is understandable. Besides, 
we note that the stress in [110] direction experiences an abrupt decrease process after strain 
up to 0.24. This is due to the fact that the corresponding four Pu— O bonds (make an 
angle of 45° with the pulling direction) have been pulled to fracture. The fracture behaviors 
have been clarified by plotting valence electron charge density maps (not shown). Under 
the same strain, the abrupt increase of spin moment can be clearly seen [Fig. EJ^c)]. While 
the spin moments in [110] and [111] directions only increase from 4.13 to 4.23 and 4.33 /i^, 
respectively, the spin moments in [001] and [100] directions are increased up to about 5.27 
and 5.11 respectively, at the end of tensile deformation. In addition, the evolutions 
of the lattice parameters with strain in Fig. [3] clearly show that along with the increase 
of the lattice parameter in the pulling direction, other two lattice parameters vertical to 
the pulling direction are decreased smoothly for tensile strains along [001], [110], and [111] 
directions. In [001] direction, the evolutions of the lattice parameters along [100] and [010] 
directions are absolutely same due to the structural symmetry. For all these three tensile 
deformations, no structural transition has been observed in our present FPCTT study. The 
structural transition over e=0A0 for [001] direction tensile deformation is mainly due to the 
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FIG. 4: (Color online) Valence charge density maps for AFM Pu02 in (Oil) plane under strains 
of (a) 0.0 and (b) 0.5 along [001] direction tensile strain. The contour lines are drawn from 0.0 to 
1.5 at 0.1 e/A 3 intervals. 

magnetic structure. 



2. Electronic properties under tension 

In order to investigate carefully the physical nature of PuC>2 under tension, in the following 
we systematically study the electronic structures of its AFM phase within LDA+ U formalism 
for pulling along the [001] direction. The valence charge density maps under strains of 
0.00 and 0.50 along [001] direction tensile deformation are plotted in Fig. HI Clearly, the 
interatomic distances in the [011] direction are elongated and that in the [100] direction are 
shortened under tension. Under strain of 0.50, fractures of the Pu— O bonds can be seen. 
These results are understandable based on our foregoing statements for its pressure behaviors 
of total energy, stress, and lattice parameters. For explicitly indicating the ionic/covalent 
character of PuC>2 under tension, we further plot in Fig. |5] the evolutions of Pu— O bond 
length in (011) plane, correlated minimum values of charge density along the Pu— O bond, 
and number of electrons transfer from each Pu to O atom under tensile deformation along the 
001 ] d irection. The electron transfer analysis is performed according to the Bader analysis 
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24J ] and similar ionic/covalent character investigations for thorium dioxide and hydrides 



have been conducted in our previous works [25|, [26[. Obviously, as indicated by Fig. E^a), 



the Pu— O bond length in (011) plane is elongated by tensile stress in [001] direction. The 
initial minimum value of charge density (0.53 e/A 3 ) along the Pu— O bond, comparable to 
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FIG. 5: (Color online) Dependence of (a) bond length of Pu— O bonds in (011) plane, (b) correlated 
minimum values of charge density along the bonds, and (c)number of electrons transfer from each 
Pu to O atom on the tensile strain along the [001] direction. 



that along the Np— O bonds included in our previous study of Np0 2 [27] . is prominently 
larger than that along the Th— O bonds (0.45 e/A 3 ) in TI1O2 25|- This indicates that the 



Pu— O and Np— O bonds have stronger covalency than the Th— O bonds. Under tension, 
this minimum value of charge density for Pu— O bond decreases near linearly to 0.26 e/A 3 at 
the end of tensile deformation, which explicitly illustrates a decreasing behavior of covalent 
feature for Pu— O bond. On the other hand, the number of electrons transfer from each 
Pu to O atom under tensile deformation also decrease with strain. Electrons of each atom 
are localized to their dominated region and thus reduces the ionicity for Pu02- Overall 
speaking, the strong mixed ionic/covalent character of Pu— O bond is weakened by tensile 
strain. 

Figure [6(a) shows the dependence of the insulating band gap E g on the tensile strain 
along the [001] direction and [6(b) presents the total DOS as well as the projected DOS 
for the Pu 5/ and O 2p orbitals under strains up to 0.00, 0.32, and 0.50. Interestingly, 
the band gap varies smoothly under strain from 0.00 to 0.12 and then becomes to decrease 
near linearly from 1.66 eV to zero under strain of 0.12 to 0.32. This explicitly indicates that 
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FIG. 6: (a) Dependence of the insulating band gap on the tensile strain along the [001] direction, 
(b) Total DOS and PDOS under strains up to 0.00, 0.32, and 0.50. The Fermi energy level is set 
at zero. 

along the [001] direction tensile deformation PuC>2 will occur an insulator-to-metal transition 
after tensile stress exceeds about 79 GPa. Figure [6](b) further illustrates this conclusion. 
The DOS at strain of 0.50 shows that the electrons of Pu and O atoms are separated from 
hybridization feature to exhibit one-atom character. And the main contribution for the 
metallic property is from both Pu 5/ and O 2p orbitals. 



IV. CONCLUSION 



In conclusion, we have studied the structural, mechanical, magnetic, and electronic prop- 
erties of PuC>2 under tension by first-principles DFT within the LDA+ U formalism. Our 
calculated ideal tensile strengths are 81.2, 80.5, 28.3, and 16.8 GPa, respectively, for pulling 
in the [001], [100], [110], and [111] directions. The [001] and [100] directions are prominently 
stronger than other two directions since that more Pu— O bonds participate in the pulling 
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process. While the spin moments in [110] and [111] directions only increase a little, the 
spin moments in [001] and [100] directions increase up to about 5.3 and 5.1 /ib, respectively. 
Through charge and density of states analysis along the [001] direction, we conclude that the 
initial strong mixed ionic/covalent character of Pu— O bond is weakened by tensile strain 
and PuC>2 will occur an insulator-to-metal transition after tensile stress exceeds about 79 
GPa. The main contribution for the metallic property is from both Pu 5/ and O 2p orbitals 
at high strain domain. 
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